arXiv: 1509.02532v3 [astro-ph.CO] 6 Apr 2016 


Prepared for submission to JCAP 


Revisiting the NVSS number count 
dipole 


Prabhakar Tiwari“ and Adi Nusser “ 

“Physics Department and the Asher Space Research Institute - Technion, Haifa 32000, Israel 
E-mail: ptiwari@physics.technion.ac.il, adi@physics.technion.ac.il 

Abstract. We present a realistic modeling of the dipole component of the projected sky distribution 
of NVSS radio galaxies. The modeling relies on mock catalogs generated within the context of ACDM 
cosmology, in the linear regime of structure formation. After removing the contribution from the 
solar motion, the mocks show that the remaining observed signal is mostly (70%) due to structures 
within z < 0.1. The amplitude of the model signal depends on the bias factor b of the NVSS mock 
galaxies. For sources with flux density, S' > 15 mJy, the bias recipe inferred from higher order 
moments is consistent with the observed dipole signal at 2.12a. Flux thresholds above 20 mJy yield 
a disagreement close to the 3a level. A constant high bias, 6 = 3 is needed to mitigate the tension to 
the ^ 2.3cr level. 
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1 Introduction 

The standard cosmological paradigm is based on the “Einstein’s Cosmological Principle” [1, 2], i.e. 
the approach to homogeneity and isotropy on the largest observable scales. Galaxies in the NRAO 
VLA Sky Survey (NVSS) [3] have been assumed to lie at redshift z 1, corresponding to a comoving 
distance of ~ 3.3Gpc for the ACDM cosmological background [4]. The galaxy clustering at these 
distances should contribute very little to sky maps of NVSS galaxy and flux distributions on large 
angular scales. The observed dipole component of these maps should be almost entirely due to Doppler 
and aberration effects [5] caused by the local solar motion as measured by the Cosmic Microwave 
Background Radiation (CMBR) [6-11]^. There have been several attempts to measure the dipole 
in radio catalogs [12-18] and somehow the dipole measured from the NVSS is puzzling with several 
independent analyses [13-18] yielding values which disagree with the CMBR velocity dipole [14, 17, 18] 
at ^ 3tT significance. 

The first extraction of the solar motion from the NVSS dipole was done in [13]. The authors 
[13] reported a velocity which was roughly in the same direction as CMBR but 1.5 to 2 times larger 
in amplitude. However, nine years later, using the same data, a huge solar motion, 1600± 400 km 
s“^, was reported [14] from both the projected number density (per steradian) and brightness dipole 
measurements. This significant (> 3cr) disagreement with the observed solar motion, immediately 
attracted the attention of the community resulting in further detailed analyses of the NVSS dipole 
[15];[16]([19]);[17, 18]. 

The authors of [15] derived a dipole almost six times larger than the CMBR velocity dipole, which 
was attributed to observational biases. Another analysis by [16] yielded a dipole with an amplitude 
roughly three times larger than the prediction from the CMBR value. A more recent analysis by [17] 
found a 2.7 — 3.2tT disagreement between the NVSS number dipole and the CMBR for various flux 
density cuts. Conclusively, the present status of NVSS dipole is a disagreement with CMBR predicted 
velocity dipole at 2.7a or more. 

However, the contribution of the intrinsic variation of the projected number density as a result 
of the NVSS galaxy clustering remains to be properly assessed. The intrinsic dipole mainly depends 
on the following key ingredients: 

^The measured CMBR dipole magnitude is of order ~ 10“® [6—11] and the corresponding speed is found to be 369± 
0.9 km s“l in the direction, I = 263.99° ± 0.14°, b = 48.26 ± 0.03° in galactic coordinates [10, 11]. In J2000 equatorial 
coordinates, the direction parameters are RA = 167.9°, DEC = —6.93°. 
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• number density distribution N(z) as a function of redshift. This dictates the fractional contribu¬ 
tion of the large scale structure from different redshifts. Available constraints on N{z) strongly 
indicate that it peaks at redshifts lower than z ^ 1. 

• relation between the clusteriirg of radio galaxies versus the underlying mass as a function of 
redshift. This biasing relation can be parametrized in terms of the galaxy, Pg{k), and matter 
power spectrum, P{k), as, Pg{k,z) = b^{z)P{k, z). The bias factor, b, could be as large as 
2 — 3 and hence produce a non-negligible boost to the dipole signal, especially if N{z) peaks at 
redshifts lower than z ^ 1. 

• constraints from the Local Universe. Our cosmological neighborhood within 100 Mpc is moving 
with ^300 km s“^ [20] and roughly aligned with the solar motion. This motion is caused by 
even farther structures, producing an alignment between the solar motion and intrinsic dipole. 

We aim at a realistic assessment of the NVSS dipole within the context of the ACDM model taking 
into account the above ingredients. In Ref. [21] (hereafter NT15), the authors approximate the 
observed redshift distribution, N°^^{z) from the best available redshift catalogs Combined EIS-NVSS 
Survey of Radio Sources (CENSORS) [22, 23] and Hercules [24, 25]. They extract the biasing factor 
b{z) by matching the ACDM predictions to the NVSS measured angular power spectrum. The redshift 
dependent biasing and number distribution are also found to be consistent with the observed radio 
luminosity functions [23, 26, 27]. The results obtained in NT15 are based on high order multipoles 
Z > 4, where the analysis is fairly insensitive to complications related to the partial sky coverage. In 
this paper we focus on the dipole, Z = I, and perform the comparison via ACDM mock catalogs. 

The outline for the paper is as following. In §2, we briefly describe the NVSS catalog. In §3.1 we 
review the basic formalism for describing the dipole from a projected map of galaxies. We provide our 
treatment for the number density N{z) and galaxy biasing in §3.2 & §3.3, respectively. §4 describes 
the NVSS mock catalogs. We present our results in §5 and conclude with a general discussion in §6. 

2 Data 

The full NVSS catalog contains ^1.7 million sources observed at I.4GHz with an integrated flux 
density S'i. 4 ghz > 2.5 mJy. The full width at half maximum angular resolution is 45 arcsec and 
nearly all observations are at uniform sensitivity. Most of the sources (90%) are unresolved. The 
catalog covers the sky north of declination —40° (J2000), which is almost 82% of the celestial sphere. 
However, the effective sky coverage for our study here is 75% since we impose the following two 
cuts: 1) we mask out sources at latitudes jfcj < 5° to avoid Galactic contamination and 2) we remove 
additional 22 sites around bright local extended radio galaxies identified by [13], to avoid overwhelming 
correlations which may leak to the large scales of interest. These are the only cuts we impose on the 
data. The authors of [17] have studied the robustness of the measured NVSS dipole with respect to 
additional cuts. They have demonstrated that masking out latitudes of < 10°, around the Galactic 
and Super-Galactic planes changes the dipole by less that 10%. The same authors have also assessed 
the effect of removing sources within 30 arcsec of known nearby galaxies listed in [28-33]. This 
procedure also resulted in small changes in the NVSS dipole. 

As shown by Ref. [13], the NVSS suffers from significant systematic mis-calibration across the 
sky at flux densities S' 15 mJy, which affects lower multipoles including the dipole. Therefore, we 
restrict our analysis to sources brighter than 15 mJy and consistently give the results for S > 15, 20, 
30, 40 and 50 mJy. Furthermore, for completeness and comparison to other authors [13, 16, 17], we 
give results for S > 10 mJy, for one set of mocks with most conservative parameters. We impose a 
maximum flux density cut S < 1000 mJy to remove extra bright sources. 

3 Basics 

3.1 The dipole component 

Let J\f{r) be the projected number density (per steradian) in the line-of-sight direction f, and M be 
the mean number density averaged over the sky. Further, let the density contrast A(f) = Af{r)/M—1. 
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Given the measured angular positions of galaxies, we employ the HEALPix [34] scheme to obtain Af{r) 
in equal area pixels. We expand A(f) in spherical harmonics, Yim(r), as 

OO ~\-l 

= X] X! ^lrnYlm{r) ■ (3.1) 

I—I rn— — l 

For a full sky coverage, the coefficients, aim can easily be calculated by inverting equation (3.1), as 

aim= f dnA{r)Yim{r). (3.2) 

J 47r 

The angular power, G/, for multipole I is, 


Cl = (jaimP) • 


(3.3) 


For a partial sky coverage, we resort to an approximate scheme [35] to compute G/, 

^obs _ _ L 

' Jim M 


(3.4) 


where is obtained from equation (3.2) but with angular integration over survey region only. 
Further, Jim = f lYimPdfl is the approximated correction suggested in reference [35]. The term 

survey 

^ describes the contribution of the Poissonian shot-noise. The G°'^® is a measure of anisotropy at 
angular scale tt jl. 


3.1.1 The solar motion contribution to the dipole 

A moving observer sees a brighter sky in the forward direction due to Doppler boosting and aberration 
effects [5]. The flux density (energy per area per frequency) of sources in NVSS data follows a power- 
law dependence on frequency S' oc v~°‘, with a « 0.75 [5]. Let v be speed of the observer relative 
to a source emitting radiation intrinsic frequency r'rest- To lowest order in u/c, the observer sees the 
radiation at a frequency z/obs = Vre&td, where 


(5 « 1 -I- (u/c) COS 0 , 


(3.5) 


and 9 is the angle between v and the line of sight to the source. Therefore, the observed flux density 
is [5], 

Sobs = Srest^^+“ (3.6) 

at a fixed frequency in observer’s frame. The extragalactic radio source population count above some 
flux density cut, follows a power law n{> S) oc S““ [5], where a; ~ 1 [17]. Therefore, the number count 
above some Sobs will change by a factor due to Doppler effect. In addition, the aberration of 

light changes the solid angle as, 

dflobs = dQ^estS~^. (3.7) 


Using these equations, the dipole in the projected number density due to the local motion, at leading 
order in v/c, is [5, 17], 

D = [2 + xil + a)]iv/c) . (3.8) 

The dipole term of the power spectrum, Gi, corresponds to to a dipole amplitude, D (|D|) as in 
equation (3.8), given by [15], 
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3.1.2 Contribution from large scale structure 

In linear theory for gravitational instability, the theoretical angular power spectrum Ci is given in 
terms of the 3D mass density power spectrum P{k) at z = 0 as [21] 

Cl = |u/m| ^ 

= ^ J dkPP{k)W^{k) , (3.10) 

where W{k) = D{z)b{z)p{r)drji{kr) is the window function in /c-space and P{k) is the ACDM 
mass density power spectrum. Further, D{z) is the linear growth factor normalized to unity at z = 0 
and p(r)dr oc 7V(z)dz is the probability for finding a galaxy in the range r — (r + dr) and r is the 
comoving distance corresponding to redshift z. 


3.2 The redshift distribution 

We rely on the results of NT 15 for biasing and the redshift distribution for radio galaxies. These 
authors modeled the redshift distribution of NVSS galaxies on the basis of the observed distribution 
iV°'^®(z) obtained from CENSORS and Hercules surveys. The CENSORS catalog was developed to 
study the evolution of the steep-spectrum radio luminosity function and the survey is presumably 
complete for flux densities S > 7.2mJy at 1.4 GHz, containing 135 sources. Indeed, the catalog covers 
only 6 deg^ of the ESO Imaging Survey Patch D (EISD). The redshift measurements of the sample 
are 73% spectroscopically complete and for the remaining sources the redshift is estimated by Kz or 
Iz magnitude-redshift relation. The Hercules redshift survey of 64 sources within 1.2 deg^ is complete 
for S > 2mJy at 1.4 GHz. We combine GENSORS-bHercules data to model NVSS source redshift 
distribution. We find 133 sources above 10 mJy flux density S. For higher flux density cuts we use 
the same as the number of sources reduces significantly and effectively there is no change in fit. We 
add that the CENSORS-bHercules combined data sets have been used to study the redshift dependent 
properties [2 1-23]. NT15 use the parametric model 


^model ^ 



(3.11) 


where the parameters oi, 02 and 03 have been fixed by fitting the CENSORS-I-Hercules redshift 
distribution as well as the NVSS angular power spectrum for I > 4. The best fit values were found to 
be oi = 0.74 ±0.57, 02 = 0.71 ±0.79 and 03 = 1.06 ±0.53. A histogram of the observed distribution, 
N°'°^{z), is plotted in figure 1 together with the model fit in equation (3.11) for the best fit values of at. 
Further we have plotted 5 curves representing V(z), for at randomly drawn from the full probability 
distribution function computed in NT15. 


3.3 Galaxy biasing 

To model the biasing of radio galaxies, NT15 use an elaborate model based on halo biasing [36] and 
the relationship between galaxy stellar masses and the prevalence of radio activities [37]. They provide 
the following convenient parametric fit. 

b^^{z) = biz^ + b2Z + bn , (3.12) 

where the best fit values of bi are bi = 0.33, 62 = 0.85 and 63 = 1.6. Furthermore, we shall employ 
other choices for the bias factor to fit with NVSS dipole observation. Particularly, we consider constant 
bias factors b{z) = 2 and 3. The thick red line in figure 2 is b^^ (z) for the best fit parameters bi, while 
the thin black lines represent random bi drawn from the full probability function from NT15. Several 
estimates of the bias factor available in the literature are also plotted in the figure. The values from 
Allison et al. [38] and Lindsay et al. [39] have a large uncertainty and agree with the bias factors 
used in this work. However, the bias of faint galaxies (Si .4 > 90 A^Jy) [40] is apparently low. The bias 
measurements in [38-40] correspond to large redshift bins and the redshift for each data points is the 
median or effective redshift of the bin. The bias of NVSS galaxies obtained by Blake et al. [41] (grey 
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Figure 1. The histogram is the observed number of sources from combined CENSORS+Hercules data sets 
with flux density S' > 10 mjy and per redshift bin = 0.2. Best fit model (input model) follows equation 
(3.11) with parameters oi = 0.74, 02 = 0.71 and as = 1.06. The number density from mock catalogs for same 
area as observation is also shown. Note that the mocks include the ACDM density fluctuation, biasing b(z) 
and the Poisson shot-noise. 


band), is assumed to be constant with redshift in their analysis and does not take into account the 
dipole. Unfortunately, we do not have any radio galaxy bias observation at redshift z 0.1 where the 
dipole signal is mostly formed. 

4 The Mock catalogs 

The origin of the intrinsic dipole signature is fluctuations on large scales with negligible coupling to 
nonlinear evolution on small scales. Hence it suffices for our purposes to construct mock NVSS catalogs 
based on linear theory for structure formation. The mocks are extracted from a random Gaussian 
realization of linear density and velocity fields sampled on a 512^ cubic grid of 17.408 comoving Gpc on 
the side. The box is equivalent, in volume, to a sphere extending to redshift z ^ 7. The grid spacing 
of 34 Mpc is well below the physical scales relevant for the dipole. The realization of the linear dark 
matter density, Sm, in the box is generated with the GRAFIC-2 [42] package and the corresponding 
velocity field is computed using the linear theory relation between velocity and density [43]. We 
adopt the latest Planck results [4] for the cosmological parameters: Hubble constant Hq =67.8 km 
s“^ Mpc“^, total matter density parameter flm = 0.308, baryonic density parameter Hb = 0.0486, 
linear clustering amplitude on 8 h~^Mpc scale, erg = 0.815 and a spectral index = 0.9667. We 
use the fitting formulae from Ref. [44] for the AGDM power spectrum. “Observers” are placed at the 
centers of spheres of 100 Mpc in radius and having with bulk motions of ~ 300 km s“^. This is 
designed to match the observed bulk flow as measured in [20] . The mean number of the galaxies in a 
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Figure 2. The redshift dependent b^j.{z) in equation (3.12) is plotted as the thick red curve for the best 
fit parameters. Thin curves correspond to random choices of the parameters appearing in bj^^(z). For 
comparison, results from the literature are also shown, as indicated in the hgure. The abbreviation FGM 
stands for ‘FIRST-GAMA-matched’. 


bin at redshift ^bin, iV(bin), is, 

A^(bin) = h(zbin) X dv X [1 + &(zbin)£>(-Zbin)^m], (4.1) 

where h(zbin) is the mean number density at the bin and dv is the bin volume. The number density 
nizhin) is obtained from the redshift distribution form the equation (3.11). The actual number of 
galaxies in the bin is a random integer, drawn from a Poisson distribution with mean 7V(bin). 

In order to include scatter due to uncertainties in the fits of N{z) and for the biasing form b^^{z) 
in equation (3.12), we draw random realizations of the parameters and bi (see equations (3.11) & 
(3.12)) from the full probability distribution for the observed {I > 4) computed in NT15 (see 
figures 1 & 2). 10^ NVSS mocks are generated, respectively, from these random sets rather than the 
mean values. The biasing form b^r^ is unsuitable for high redshifts but this is insignificant since the 
contribution to the dipole from large scale structure at z > 1 practically vanishes (see figure 3). For 
computational purposes we simply take b{z > 1.5) = b{z = 1.5). When exploring constant bias values, 
we still generate 10^ mocks for random shot noise and Oi. The angular orientation of the observed 
NVSS sky in relation to the direction of the observed dipole is matched exactly. We rotate the mocked 
map such as the dipole direction matches with the observed direction and then we mask the Galactic 
plane and 22 bright local extended sites as identified in [13]. These final mock NVSS galaxy maps are 
used to calculate Ci and are compared with NVSS observed dipole. 

It is important to explore the range of redshifts and scales of structures which are mostly re¬ 
sponsible for the intrinsic dipole signature. In figure 3 we show the fractional contribution to Ci from 


- 6 - 















structure up to a certain redshift from the theoretical relation (see equation (3.10)). The dipole is 
mostly determined by low redshift structure with 60% of the power already achieved by z = 0.1. 

In hgure 4 we show the fc-space window function, W{k)^ and the theoretical differential contri¬ 
bution Cl to per log(fc) bins. Although the window function peaks around k ^ 0.001 hMpc“^, the 
overall dipole signal gets most of its contribution from scales k ~ 0.01 hMpc“^. 



Figure 3. The fractional contribution to Ci from structures below redshift z. 


5 Results 

The statistical assessment of the ability of the models to match the observed is done on the basis 
of counting the number of mocks with Ci exceeding the observed value. We perform this analysis 
for several lower cuts, 5'min, on the flux density, S, where for each cut we use 10^ random mocks. 
The observed is calculated after removing the contribution of the solar motion from the NVSS 
number density maps, according to the relation (3.8). 

We first explore consequences of b{z) = b^^{z) in (3.12) with the parameters which yield a 
good match to Ci for I > 4 (see NT15). Figure 5 shows the distribution of Ci in the 10® mocks 
for 5'inin = 20 mjy. The distribution is well approximated by a distribution with 3 degrees of 
freedom corresponding to the average of common variance of three independent Gaussian variables 
(< (|aiip + |aioP + |aiiP >)• The theoretical distribution is also plotted in the figure 5, as the 
red curve. The probability for the mocks to yield a Ci larger than the observed NVSS dipole (7°®® is 
0.2% for b^^. This corresponds to ~ 2.89 ct significance level for a normal distribution of one variable. 
Values of for several flux thresholds, S'min, are summarised in table 1. Note that the NVSS dipole 
for S'min = 10 mJy may suffer from systematic biases [13] and we consider it unreliable. We list results 
with this cut only for the sake of completeness and comparison with other authors. According to the 
results in table 1, the ACDM mocks with the bias b^^ are reasonably consistent at the 2.12a level with 
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Figure 4. The dipole fc-space window function, W{k), and the theoretical Ci par log{k) bins, 
integral over per log{k) is normalised to ACDM Ci. The dipole (Ci) is most sensitive to fc ~ 0.01 


The ACi 
hMpc"\ 


the NVSS dipole with Smin = 15 mJy. Nonetheless, there is a tension between the mocks and the data 
at higher fluxes. The biasing template is inferred from high order multipoles which are sensitive 
to the galaxy distribution at redshifts z ^ 0.2 (see figure 7 in NT15). The dipole signal, however, is 
largely generated at lower redshifts, z < 0.1 by fluctuations with wave numbers k ^ O.OlhMpc”^ (cf. 
figures 3 & 4). Therefore, the choice b{z) = b^^ may not be suitable for the dipole in the case of a 
redshift and scale dependent (radio) galaxy biasing. Alternatively, we generate mocks for b{z) = 2 
and 3. The results for multipoles I = 1,2,3 and 4 are shown in table 2. Both values of constant b are 
consistent with the dipole for S'nUn = 15 and, as with b^^, higher flux thresholds are in tension with 
the observed dipole for 6 = 2. The high value 6 = 3 tones down the dipole tension to an acceptable 
level. However, this value is significantly high compared to most observational estimates [21, 38-41] 
as seen in figure 2 . 


^min (mJy) 

(RA, Dec) 

il,b) 

(xlO^) 

shot-noise (xlO®) 

p-value 

a- significance 

10 

(135°, 23°) 

(204°, 38°) 

0.82 

1.8 

0.02389 

1.98 

15 

(151°,-6°) 

(246°, 38°) 

1.10 

2.5 

0.01696 

2.12 

20 

(151°, -14°) 

(253°, 32°) 

1.96 

3.3 

0.00193 

2.89 

30 

(160°,-16°) 

(263°, 36°) 

2.37 

4.7 

0.00262 

2.79 

40 

(150°, -35°) 

(268°, 16°) 

2.89 

6.3 

0.00173 

2.92 

50 

(174°, -36°) 

(286°, 24°) 

3.79 

8.0 

0.00209 

2.86 


Table 1. The observed and shot-noise (^) along with dipole direction. The probability is the count of 
mocks with Ci exceeding the observed The corresponding a- significance is also given for the observed 

probability. The p-value is the probability (Ci > (7°'^'*). 












Smin (mJy) 

Cl 

C 2 

C 3 

C 4 


p-value 

cr-significance 

p-value 

cr-significance 

p-value 

(T-significance 

p-value 

(T-significance 





q 

II 

-0 





15 

0.01437 

2.19 

0.00324 

2.72 

0.32412 

0.46 

0.12362 

1.16 

20 

0.00149 

2.97 

0.00493 

2.58 

0.96881 

— 

0.03533 

1.81 

30 

0.00300 

2.75 

0.00765 

2.43 

0.69182 

— 

0.30641 

0.51 

40 

0.00252 

2.80 

0.00696 

2.46 

0.41254 

0.22 

0.41076 

0.23 

50 

0.00154 

2.96 

0.08581 

1.37 

0.46751 

0.08 

0.31316 

0.49 





b{z) = 3.0 





15 

0.06170 

1.54 

0.02345 

1.99 

0.56519 

— 

0.34328 

0.40 

20 

0.00922 

2.36 

0.02290 

2.00 

0.98633 

— 

0.13678 

1.10 

30 

0.01011 

2.32 

0.02437 

1.97 

0.80103 

— 

0.49466 

0.01 

40 

0.00656 

2.48 

0.01882 

2.08 

0.53594 

— 

0.56074 

— 

50 

0.00391 

2.66 

0.13536 

1.10 

0.57366 

— 

0.43821 

0.16 


Table 2. The low multipole Ci, C 2 , C 3 and Ci for different biasing templates. The p-value is the probability 
[Cl > For all p-values (< 0.5) we have given the corresponding cr-significance for a normal distribution 

of one variable. 



C, 


Figure 5. The mocks Ci distribution and the observation probability from 1,00,000 samples is shown for 
the flux density cut S' > 20 mjy. The dotted line is mean Ci. Note that the y-axis is in log scale as the 
probability for NVSS observed dipole or a dipole with larger amplitude is extremely tiny (0.2%) and not even 
visible in linear scale. A good fit to mocks Ci distribution to X^(3) probability distribution is also shown. 
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6 Conclusion and Discussion 


We have presented an analysis of the NVSS dipole incorporating, for the first time, the contribution 
from the large scale inhomogeneities in the galaxy distribution. Mock catalogs mimicking the NVSS 
catalog are used to model uncertainties in the galaxy redshift distribution and incomplete angular 
coverage of the data. The bias factor obtained from / > 4 by NT15 is at odds with the number 
count dipole measurement at the ^ 2.9cr significance level for the ACDM cosmology. However, the 
intrinsic dipole signal is shapes at low (z ^ 0.1) redshifts and on large scales (k ~ O.OlhMpc"^); about 
70% (45%) of the dipole (I = 4) contribution is from redshift 2 < 0.1. Therefore, it is prudent to 
consider a large bias factor for radio galaxies at low redshifts as an explanation to the dipole signal. 
Independent bias estimates at such low redshift are crude, as the redshift of radio galaxies is generally 
unknown and we do not have much galaxies in low redshift bins. A direct assessment of the redshift 
dependence of the bias of radio galaxies is done in wide redshift bins due to the small number of radio 
galaxies with redshift information. A summary of some bias factor determination in the literature is 
given figure 2. The radio galaxy bias at low redshift z < 0.1 is loosely constrained and if the bias 
is as high as b(z < 0.1) « 3, it tones down the disagreement between the model and the data to a 
^ 2.3cr significance level for S > 20, 30, 40 and 50 mJy. This bias value agrees, within the large error 
bars, with the estimate obtained in Ref. [39] (see figure 2). But it is substantially high with respect to 
other results [38, 40, 41], including the value obtained by NT15 for the same NVSS data from higher 
multipoles I > 4. We emphasize that the bias extracted in NT 15 follows the detailed halo biasing 
recipe [36] and a relation between galaxy stellar mass and radio activity [37]. 

The NVSS galaxy sky distribution may be affected by the presence of a large scale mass dis¬ 
tribution which generates the observed bulk motion of a sphere of radius 100ft.“^ Mpc around the 
observer. The modeling here incorporates the presence of this mass fluctuation by selecting “ob¬ 
servers” satisfying the bulk motion constraint. We also rotate the maps so that the dipole signal is 
aligned with the observation mask as in the real data. As seen in table 3, these two constraints give 
modest improvement in matching the model to the observations. 

There are certain claims in the literature for excess of power on scales of 100s of Mpcs relative to 
the ACDM model [45-47]. We will refrain from a detailed assessment of these claims and we simply 
emphasize the following point. The dipole signal is matched for the ACDM for 6 = 3. Therefore, any 
such power excess must be 3 fold higher than the ACDM. Such large excess at 2 ~ 0.1 will generate 
a bulk flow within a 1006.“^ Mpc which is 3 times the predicted value of the bulk flow in the ACDM, 
i.e. a motion of ~ 800 kms“^. Such a value of the bulk flow is at odds with the highest bulk flow 
estimates in the literature. 

There are several non-conventional explanation for the NVSS dipole. For example, the anisotropic 
mode generated before inflation may re-enter the horizon at late time and produce large scale anisotropy 
[48-50]. Alternatively, the primordial non-Gaussianities (NG) may affect the structures at large scales. 
Particularly, significant NG in inflationary cosmology scenario may potentially give an increase in bi¬ 
asing over large scale and produce excess dipole signal [51]. In Ref. [52] the authors argue that 
NVSS dipole anisotropy can be partially explained with a local over- or under density volume. The 
authors in Ref. [52] constrain themselves to local structures consistent with GMBR dipole and other 
observations. 

The NVSS dipole remains puzzling within the ambit of radio galaxy biasing. The extra-large 
power required to explain dipole signal is ^ 3 fold in comparison to the AGDM and very unlikely. 
The NVSS excess dipole is most likely due to a residual calibration error across the sky which also 
affects the high flux cuts. The radio galaxy dipole observation will immensely improve with the Square 
Kilometer Array. 
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